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We consider the "Touch and Stop" cluster growth percolation (CGP) model on the two dimen- 
sional square lattice. A key-parameter in the model is the fraction p of occupied "seed" sites that 
act as nucleation centers from which a particular cluster growth procedure is started. Here, we 
consider two growth-styles: rhombic and disk-shaped cluster growth. For intermediate values of p 
the final state, attained by the growth procedure, exhibits a cluster of occupied sites that spans the 
entire lattice. Using numerical simulations we investigate the percolation probability and the order 
parameter and perform a finite-size scaling analysis for lattices of side length up to L = 1024 in 
order to carefully determine the critical exponents that govern the respective transition. In contrast 
to previous studies, reported in [Tsakiris et at, Phys. Rev. E 82 (2010) 041108], we find strong 
numerical evidence that the CGP model is in the standard percolation universality class. 
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I. INTRODUCTION 

The pivotal question in standard percolation [l], Q is 
that of connectivity. A basic example is 2D random 
site percolation, where one studies a lattice in which a 
random fraction p of the sites is "occupied". Clusters 
composed of adjacent occupied sites are then analyzed 
regarding their geometric properties. Depending on the 
fraction p of occupied sites, the geometric properties of 
the clusters change, leading from a phase with rather 
small and disconnected clusters to a phase, where there 
is basically one large cluster covering the lattice. Therein, 
the appearance of an infinite, i.e. percolating, cluster is 
described by a second-order phase transition. 

There is a wealth of literature on a multitude of vari- 
ants on the above basic percolation problem that model 
all kinds of phenomena, ranging from simple configura- 
tional statistics to "string" -bearing models that also in- 
volve a high degree of optimization, e.g. describing vor- 
tices in high T c superconductivity [3|, |4J] and domain wall 
excitations in disordered media such as 2D spin glasses 

@ and the 2D solid-on-solid model Q ■ Besides dis- 
crete lattice models there is also interest in studying con- 
tinuum percolation models, where recent studies reported 
on highly precise estimates of critical properties for spa- 
tially extended, randomly oriented and possibly overlap- 
ping objects with various shapes Q. 

One such variant of the above basic percolation model 
is the recently proposed "Touch and Stop" cluster growth 
percolation (CGP) model 0,113] • In the CGP model, a 
random fraction of p sites is distinguished to comprise a 
set of "seed" sites for which a particular cluster growth 



'Electronic address: olivcr.mclchert@uni-oldenburg.dc 



procedure is evolved (see sect.fll]). I.e., starting from the 
seed sites, clusters are grown by assimilating all near- 
est neighbor sites layer by layer in an iterative, discrete 
time fashion. As soon as one cluster comes in contact 
with other clusters, the growth procedure for all involved 
clusters is stopped. When there is no growing cluster 
left, the cluster growth procedure is completed and the 
connected regions of adjacent occupied sites, i.e. the final 
clusters, are analyzed regarding their geometric proper- 
ties. To support intuition: if p is rather small, the growth 
of a particular cluster is unlikely to be hindered by other 
clusters within the first few time steps, since the typ- 
ical distance between the respective seed sites is large 
compared to the lattice spacing. Consequently, the final 
clusters are rather large, see Fig. HJa). As p increases, 
the typical size of clusters in the final configuration de- 
creases due to the increasing density of initial seed sites, 
see Figs. QJa-c). Increasing p even further leads to an 
increasing size of the largest cluster. This is due to the 
larger probability to yield adjacent seed sites already in 
the initial configuration, for which a growth procedure 
is subsequently inhibited. For increasing p and due to 
spacial homogeneity this in turn leads to an increase of 
the size of the largest cluster, see Figs. Hfc— e). 

The two distinct scaling regimes of the largest cluster 
size were discussed in Ref. @ . The behavior at small val- 
ues of p was illustrated to be a finite size effect and the 
scaling at intermediate values of p was found to signal a 
continuous phase transition qualitatively similar to the 
standard percolation transition. In Ref. the critical 
properties of the latter transition were addressed. Albeit 
the authors of Ref. [l(| perform a scaling analysis to es- 
timate the critical point and a set of critical exponents 
that govern the phase transition, we are convinced that 
the analysis can be improved in various ways. In this 
regard, we here revisit the continuous phase transition 



2 





(c) 




(d) 



Mm 




FIG. 1: Samples of final configurations obtained for the 
CGP model on 2D square lattices of side length L = 128. 
The snapshots relate to different values of the initial density 
of active seed sites: p = 0.001, 0.01, 0.02, 0.05, 0.495, 0.7, for 
subfigures (a-f ) , respectively. In the limit of large system sizes 
and above the critical point p c m 0.495, the largest cluster in 
the final configuration spans the lattice along at least one 
direction. In each subfigure, the largest cluster is colored 
black and all finite clusters are colored grey. 



found for the CGP model at intermediate values of p and 
perform a finite size scaling analysis considering geomet- 
ric properties of the largest cluster for the individual final 
configurations obtained using computer simulations. To 
this end, the results we report in sect. lllll below (obtained 
using two different cluster growth styles, i.e. rhombic and 
disk-shaped cluster growth) are different from those pre- 
sented in Ref. [lOj in the following respect: in contrast to 
their findings we here present strong numerical evidence 
that the percolation transition of the CGP model at in- 
termediate values of p is in the same universality class 
as standard percolation. By employing the data collapse 
technique for different observables we aim to yield maxi- 
mally justifiable results due to a high degree of numerical 
redundancy. Further, we analyze the distribution of clus- 
ter sizes right at the critical point in order to numerically 
assess a particular critical exponent that was previously 
only computed via scaling relations. 

The remainder of the presented article is organized as 
follows. In section [TTl we introduce the model in more 
detail. In section IIII1 we outline the data collapse tech- 
nique and we list the results of our numerical simula- 
tions. Therein, the discussion is focused on rhombic clus- 
ter growth. Section IIVI concludes with a summary. Ap- 
pendix [5] shows a further analysis of the order parameter 
resembling the one presented in Ref. [lOj, and appendix IB1 
elaborates on the results obtained for disk-shaped cluster 
growth. 



II. MODEL AND ALGORITHM 



Here, we consider 2D square lattices of side length L 
and N — L 2 sites. Initially, a starting configuration is 
prepared by occupying a random fraction p of the lat- 



tice sites via seed sites. These seeds indicate the cen- 
ters of a set of active clusters that will evolve during 
the growth procedure. The remaining sites are consid- 
ered empty. For the set of active clusters, the growth 
procedure consists of the repeated application of the fol- 
lowing two steps: (i) to perform a single sweep, consider 
the still active clusters sequentially in random order. For 
each cluster perform the cluster update below, (ii) so as 
to perform a single cluster update, consider the surface of 
the respective cluster, i.e. the set of nearest neighbor sites 
of those sites that build up the current cluster and that 
do not belong to the cluster, yet. To complete a cluster 
update, the "Touch and Stop" CGP model discriminates 
the following two situations: (ii. 1 ) if alien sites, i.e. sites 
belonging to clusters different from the current one are 
found among the surface sites, the cluster growth pro- 
cedure for all involved clusters is stopped. The respec- 
tive clusters are further deleted from the set of active 
clusters; (ii.2) if all surface sites are empty, amend the 
current cluster by the set of surface sites. The above 
two steps (i) and (ii) are iterated until no active cluster 
is left and the final configuration is reached. I.e., the 
evolution of the system proceeds in discrete time-steps, 
where within one particular time-step each active cluster 
is updated once. Note that the precise configuration of 
occupied sites in the final configuration might depend on 
the order in which active clusters are picked in step (i). 

Fig. [5] illustrates the cluster growth procedure for a 
square lattice of side length L = 12, where initially 18 
seed sites are present (see Fig. (21(a)), i.e. p = 0.125. The 
18 seed sites are arranged in 10 clusters. Consequently, 
the initial configuration exhibits small clusters of adja- 
cent occupied sites for which, according to step (ii) above, 
not a single cluster update is performed. Both of these 
clusters occupy 5 sites and one of them is picked to sig- 
nify the largest cluster on the lattice. In the figure it 
is colored black to distinguish it from the other clusters 
(which are colored grey) . While sweeping over the active 
clusters in step (i), a successful cluster update for one 
cluster might hinder the growth of another cluster that 
still has to be considered in the respective sweep (see, 
e.g., the two occupied next nearest neighbor sites in the 
upper left corner of Fig. [21(a))- After one sweep the con- 
figuration of occupied clusters has changed, see Fig. [2Kb). 
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FIG. 2: Example of the cluster growth procedure for a 2D 
square lattice of side length L — 12. For the initial configu- 
ration of active seed sites shown in subfigure (a), the growth 
procedure is completed after the two steps (a-b) and (b-c), 
already. Subfigure (c) indicates the final configuration. Sec 
text for more details. 
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Now, the largest cluster comprises 15 sites and there are 
only three active clusters left. The next sweep, resulting 
in the final configuration shown in Fig. Etc) , completes 
the growth procedure. Therein, the largest cluster occu- 
pies 53 sites and it spans the lattice along the vertical 
direction since its projection on the independent lattice 
axis covers 12 sites in the vertical direction. 

As pointed out in the introduction and discussed in 
Ref. Q, the model features two transitions: a sharp tran- 
sition at a very small value of p, which is due to finite 
size effects, and a continuous transition at intermediate 
values of p. Below, we will perform simulations on square 
systems of finite size and for different values of p in order 
to assess the critical properties of the CGP model in the 
vicinity of the critical point in 2D. The observables we 
consider are derived from the set of clusters in the final 
configuration, attained after the growth process is com- 
pleted. The observables will be introduced and discussed 
in the subsequent section. 

III. RESULTS 

So as to assess the critical properties of the CGP 
model in the vicinity of the critical point in 2D we 
performed simulations for square lattices of side length 
up to L = 1024. For each system size considered, we 
recorded data sets at 24 supporting points in a p-range 
that encloses the critical point on the p-axis. Each data 
set comprises the set of clusters in the final configura- 
tion for rs 10 4 individual samples. As stressed in the 
introduction, we consider two different cluster growth 
styles: rhombic and disk-shaped cluster growth, defining 
the rhombic CGP model (CGP-R) and the disk-shaped 
CGP model (CGP-D), respectively. In the remainder of 
this section we report on the results for rhombic cluster 
growth. The results for disk-shaped cluster growth are 
detailed in appendix [B] 

Most of the observables we consider below can be 
rescaled following a common scaling assumption (formu- 
lated below for a general observable y(p, L)). This scaling 
assumption states that if the observable obeys scaling, it 
can be expressed as 

y{p,L)=L~ h f[(p- Pc )L% (1) 

wherein a and b represent dimensionless critical expo- 
nents (or ratios thereof, see below), p c signifies the criti- 
cal point, and /[•] denotes an unknown scaling function. 
Following Eq. (fT]), data curves of the observable y(p,L) 
recorded at different values of p and L collapse, i.e. fall 
on top of each other, if y(p,L)L b is plotted against the 
combined quantity e = (p — p c )L a and if the scaling pa- 
rameters p c , a and b that enter Eq. ((T|) are chosen prop- 
erly. The values of the scaling parameters that yield 
the best data collapse signify the numerical values of 
the critical exponents that govern the scaling behavior 
of the underlying observable y(p,L). In order to obtain 
a data collapse for a given set of data curves we here 




FIG. 3: Finite size scaling of the spanning probability P(p) 
for the CGP-R model on 2D square lattice of side length 
L = 192 through 1024, averaged over different initial seed 
configurations. The main plot shows the data collapse ob- 
tained according to Eq. fj] and the inset illustrates the raw 
data close to the critical point. 



perform a computer assisted scaling analysis, see Refs. 
[Til [l2j ■ The resulting numerical estimates of the criti- 
cal point and the critical exponents for the 2D "Touch 
and Stop" model implemented using rhombic and disk- 
shaped cluster growth as well as the critical properties of 
the standard 2D site percolation model (for comparison) 
are listed in Tab.fl] Below, we report on the results found 
for different observables: 

a. Spanning probability: As a first observable we 
consider the probability P(p) that the final configuration 
features a cluster that spans the lattice along at least 
one direction, averaged over different initial seed config- 
urations. It is expected to scale according to Eq. ([T]), 
where a = \jv and 6 = 0. A data collapse performed 
for the four largest system sizes L = 384, 512, 768, 1024 
in the range e £ [—2, 2] on the rescaled p-axis yields the 
scaling parameters p c = 0.494(1), and v = 1.34(7) for 
a quality S = 0.71 of the data collapse, see Fig. [3] As 
it appears, the numerical value of v agrees well with the 
value for standard 2D percolation, i.e. v vcxc — 4/3. While 
the location of the critical point is in agreement with the 



TABLE I: Critical properties that characterize the CGP phe- 
nomenon (CGP-R: rhombic clusters; CGP-D: disk-shaped 
clusters) in 2D. From left to right: Critical point p c , crit- 
ical exponents v, f3, 7, and r. For comparison, the critical 
properties for the standard site percolation (SP) problem on 
the 2D square lattice (figures up to the third decimal) are 
listed also. 





Pc 


V 


P 


7 


r 


CGP-R 
CGP-D 
SP 


0.4938(7) 
0.4978(5) 
0.593 


1.33(3) 
1.32(4) 
1.333 


0.141(3) 
0.145(4) 
0.139 


2.38(5) 
2.38(3) 
2.389 


2.057(1) 
2.050(6) 
2.055 
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value reported in Ref. [l(| , the critical exponent is rather 
different (Rcf. 10] reports v = 1.17). 

b. Order parameter: As a second observable we con- 
sider the order parameter, given by the relative size 



-Pmax(p) = (Smax(p)} 



(2) 



of the largest cluster in the final configuration, averaged 
over different initial seed configurations. This observable 
scales according to Eq. ([1]), where a = 1/v and b = —/3/v. 
Again, a data collapse for the four largest system sizes 
in the range e € [—2,2] yields the scaling parameters 
p c = 0.4937(2), v = 1.34(5), and /3 = 0.141(3) for a 
quality S = 0.77 of the data collapse, see Fig.[4|b). Note 
that the numerical value of the order parameter exponent 
(3 is in agreement with the standard percolation exponent 
/3 pcrc = 5/36 « 0.139 (Ref. pjj reports (3 = 0.24). 

A different way to analyze the same set of data is in 
terms of the Binder ratio 1131 



b(p) 



(Smaxb)) 1 



(3) 



This observable scales according to Eq. ([T]), where, as for 
the spanning probability above, a = \jv and 6 = 0. The 
best data collapse yields p c = 0.4937(4), and v = 1.34(8) 
with a quality S = 0.96 in the range e 6 [—2, 2], see Fig. 
Ha). 

A further critical exponent can be estimated from the 
scaling of the order parameter fluctuations x(p)j defined 
as 



X (p) = N[(s 2 m&x (p)) 



<(P)> S 



(4) 



The order parameter fluctuations are expected to scale 
according to Eq. ([1]), where a = and b — —j/v. 

A best data collapse is attained for p c = 0.4938(7), v = 
1.33(3), and 7 = 2.38(5) with a quality S — 0.84, again in 
the range e G [—2,2], see Fig. Htc). Here, the numerical 
value of the fluctuation exponent 7 is in agreement with 
the standard percolation exponent 7 p0 rc — 43/18 « 2.389 
(Ref. [13 reports 7 = 1.91). 

c. Average size of the finite clusters: As a third ob- 
servable we consider the average size (S& n (p)} of all finite 
clusters for a particular final configuration, averaged over 
different initial seed configurations @, [l4[ • This observ- 
able is defined as 



Sfin(p) 



El g2 n s (p) 



(5) 



where n s (p) indicates the probability mass function of 
cluster sizes for a single sample at a given value of p. 
Note that the sums run over finite clusters only 0, H3] 
(indicated by the prime), i.e. if a final configuration fea- 
tures a system spanning cluster, this cluster is excluded 
from the sums. The average size of all finite clusters is 
expected to scale according to Eq. ([1]), where a = l/u and 
b = —7/^. For this observable, considering the system 
sizes L = 128, 192, 256, 384, a best data collapse is found 
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FIG. 4: Finite size scaling analyses related to the relative 
size s max of the largest cluster for the CGP-R model on 2D 
square lattice of side length L = 192 through 1024, averaged 
over different initial seed configurations. The main plots show 
the data collapse obtained according to Eq. [T] and the insets 
illustrate the raw data close to the critical point. The sub- 
figures show different ways to analyze s ma x in terms of (a) 
the Binder parameter b(p), (b) the order parameter s ma x(p) 
itself, and, (c) the fluctuation x{v) = A r var(s max ) of the order 
parameter. 
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FIG. 5: Scaling analysis of different observables at the criti- 
cal point p c of the CGP-R model. The main plot shows the 
probability mass function n a of clusters of size s, displayed 
using logarithmic binning. The inset illustrates the scaling of 
the average size (S max ) of the largest cluster and the order 
parameter fluctuations \- 



for p c = 0.491(1), v = 1.35(2), and 7 = 2.37(3) with a 
quality S — 1.81 of the data collapse. Considering the 
sequence of larger system sizes L = 384, 512, 768, 1024, 
the optimal choices for the scaling parameters result in 
p c = 0.4930(5), v = 1.35(2), and 7 = 2.35(4) with a 
quality S — 1.16. Both results agree within errorbars 
and were obtained in the range e G [— 1, 1]. Further, the 
value of 7 estimated from the average size of the finite 
clusters support the numerical value estimated from the 
scaling behavior of the order parameter fluctuations. 

d. Further observables recorded at p c : The various 
estimates for the critical point are in agreement with 
p c = 0.494(1). At this critical point, Eq. ([]} reduces 
to y(p c , L) oc L~ b . We performed additional simula- 
tions at p c to determine the scaling dimension of the 
largest cluster S max , defining the fractal dimension df 
according to (S max ) oc L df . Considering the scaling form 
(Smax)(L) = aL d f for a fit in the interval L G [20, 1024], 
we find a = 0.669(1), and df = 1.8995(4) (where the 
reduced x 2 -value reads x 2 e d = 0-65), see inset of Fig. [5] 
Albeit the error bars obtained from the least-squares fit 
is notoriously small, the numerical value of the fractal 
dimension itself is in good agreement with the ordinary 
percolation estimate c?/, pC rc = 91/48 ~ 1.896 (Ref. [10 ] 
report df = 1.79(1)). 

For numerical redundancy we further estimate the ex- 
ponent ratio 7/1/ from the order parameter fluctuations, 
Eq. (Q}, at p c via x{L) K L 1 ^, where a fit to the interval 
L G [20, 1024] yields 7/1/ = 1.772(3) (x r 2 cd = 1.0). Allow- 
ing for a slight deviation from the pure power law form 
using the effective scaling form x(L) oc (L + AL) 1 ^ , we 
obtain 7/^ = 1.785(8) (x 2 cd = 0.81), see inset of Fig. 
[SJ Note that the resulting exponent ratio is not only in 
good agreement with the numerical values for 7 and v 



obtained above, but also with the standard percolation 
estimate 7 per c /z/perc = 129/72 « 1.792 (Ref. [lOj report 
7/1/ = 1.63(1)). 

Further, at the critical point we recorded the proba- 
bility mass function n s for clusters of size s. It exhibits 
an algebraic decay following n s oc s~ r , for which we es- 
timate the exponent r = 2.057(1). This estimate was 
obtained for the data recorded at L = 1024 using a least- 
squares fit to the interval s G [200, 10000] (the errorbar 
was estimate by bootstrap resampling of the underlying 
datasets), see main plot of Fig. [5] Note that the value 
is reasonably close to the standard percolation estimate 
Tpcrc = 187/91 « 2.055 (Ref. [nj report t = 2.08, com- 
puted using a hyperscaling relation). 

Finally, at p c and on a large lattice of L = 1024 we 
estimate the fraction of occupied sites in the final con- 
figuration to be pfi n = 0.589697(5) (for a smaller lat- 
tice setup with L = 128 we find p fin = 0.58964(4)). 
These values are only slightly smaller than the fraction 
Pc, P crc = 0.59274621(13) of occupied sites at the critical 
point of 2D site percolation (l~5l. [l6j. 



IV. CONCLUSIONS 

We have investigated the continuous transition of the 
"Touch and Stop" cluster growth percolation model at in- 
termediate values of the initial fraction p of seed sites by 
means of numerical simulations, implementing rhombic 
and disk-shaped cluster growth. Previously, the critical 
point and the critical exponents that govern the transi- 
tion were estimated [l(| and it was concluded that the 
transition is in a different universality class than stan- 
dard percolation. From a point of view of data analy- 
sis, the previously presented analysis could be improved 
in various regards. Here, we have revisited the clus- 
ter growth percolation model and performed a finite-size 
scaling analysis considering the geometric properties of 
the largest clusters for the individual final configurations 
attained after the growth procedure is completed. Us- 
ing large system sizes and appropriate sample sizes, we 
obtained highly precise estimates for the critical points 
of both cluster growth styles and the usual critical expo- 
nents that characterize the scaling behavior of the span- 
ning probability, the order parameter (and its fluctua- 
tions), and the probability mass function of cluster sizes 
right at the critical point. We further used different ob- 
servables to estimate individual exponents. E.g., in order 
to determine the fluctuation exponent 7, we considered 
the order parameter fluctuations x (using the data col- 
lapse techniques shown in Fig. [He) as well as the finite- 
size scaling right at p c illustrated in the inset of Fig. [S]) 
and the average size of the finite clusters Sfi n . This lead 
to highly redundant numerical estimates for the critical 
exponents and yields maximally justifiable results. 

In summary, we find that the critical exponents esti- 
mated for the "Touch and Stop" cluster growth model 
are in reasonable agreement with those that describe 
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the standard 2D percolation phenomenon, see Tab. |U 
Hence we conclude that the 2D cluster growth percola- 
tion model is in the same universality class as standard 
2D percolation. 

Further, we found that in the vicinity of the critical 
point, the cluster growth procedure (for L > 256) is com- 
pleted after 3 — 4 steps. I.e., in comparison to ordinary 
site percolation, the growth procedure can effectively be 
seen as a process that affects the cluster configurations 
close to p c only locally. This might render it somewhat 
intuitive that the cluster growth procedure leads to a shift 
of the (anyway non-universal) critical point but does not 
change the critical exponents as compared to standard 
percolation. 
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Appendix A: Order parameter analysis similar to 
that of Tsakiris et al. 



An admissible criticism of the analysis presented ear- 
lier in sect. IIIII is, by considering the order-parameter 
and using the data-collapse method, one has in princi- 
ple three parameters (p Cl v and f3) to adjust in order to 
yield a data collapse. From a practical point of view, 
one might first consider a dimensionless quantity as, e.g., 
the percolation probability or the binder ratio, to arrive 
at estimates for p c and v (which are still two parame- 
ters to adjust at once). These can then be inserted into 
the scaling relation for the order-parameter, leaving only 
the scaling parameter /3 left to adjust in order to ob- 
tain a best collapse of the order-parameter curves. We 
found that whatever protocol we followed, the numerical 
estimates of the scaling parameters did vary only within 
their respective errorbars. Hence, from a point of view of 
analysis using the data-collapse method, the numerical 
estimates of p c , v, and /3 (as listed in Tab. IJ) appear to 
be consistent and reliable. 

However, we also performed an analysis of the order- 
parameter similar to that of Tsakiris et al., cf. Fig. (4) of 
Ref. [l0(. The benefit of their "crossing-point" method 
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FIG. 6: Order parameter analysis for the CGP-R model, sim- 
ilar to that performed by Tsakiris et al. (a) relative size s ma x 
of the largest cluster for different system sizes scaled by a fac- 
tor of L^l" ', where f3/v is adjusted such that the data curves 
exhibit a common crossing point (see text for more details), 
(b) results of a bootstrap resampling analysis for the analy- 
sis in (a). The central figure shows the distribution of tuples 
(P/v,Pc) in the plane and the small adjacent plots show the 
probability density function along the independent directions. 
The points with error bars in the central figure show the re- 
sulting estimates for the corresponding resampled parameters. 



is that it features only one scaling parameter, namely 
the exponent ratio j3/v fl7j |. In the remainder of this 
appendix we explain the "crossing-point" analysis proce- 
dure and report on the results obtained therewith. 

Starting with the data curves for the relative size s max 
of the largest cluster for different system sizes, we scale 
the data curves by a factor of of L^' v . If wc consider 
N data curves subject to j3/v > 0, we yield M — 
(N(N — l))/2 presumably different crossing points pi, 
i = 1, . . . , M. The exponent ratio /3/u is then adjusted 
in order to minimize the width Ap = max(p;) — min(pj) 
of the set of crossing points, resulting in a tuple ((3/v,p c ) 
of "optimal" parameters (we set p c = 'YlnPi/M). In or- 



7 




FIG. 7: Samples of final configurations obtained for the 
CGP-D model on 2D square lattices of side length L = 128 us- 
ing disk-shaped cluster growth. The snapshots relate to differ- 
ent values of the initial density of active seed sites: p — 0.001, 
0.01, 0.02, 0.2, 0.495, 0.7, for subfigures (a-f), respectively. 
In the limit of large system sizes and above the critical point 
p c « 0.497, the largest cluster in the final configuration spans 
the lattice along at least one direction. In each subfigure, 
the largest cluster is colored black and all finite clusters are 
colored grey. 

der to estimate errors for these optimal parameters, we 
performed a bootstrap resampling procedure of the raw 
data and analyzed each resampled set of data curves us- 
ing the above procedure. We considered an overall num- 
ber of N — 6 data curves and m = 2000 bootstrap 
samples. Fig. EJa) shows one of these bootstrap sam- 
ples. The analysis of all m optimal tuples (/3/u,p c ) in 
the plane is summarized in Fig. [SJb). Therein, the cen- 
tral plot shows the distribution of tuples in the plane 
and the small adjacent plots show the probability den- 
sity function along the independent directions. The re- 
sulting estimates for the resampled optimal parameters 
read fi/v = 0.105(11) and p c = 0.4937(8). Note that 
these are in excellent agreement with the results obtained 
from the data collapse analysis, listed in Tab. fl] Further, 
Pearson's correlation coefficient between the two param- 
eters reads p(l3/u,p c ) = —0.988(1). This indicates that 
for an increasing numerical value of the exponent ratio, 
the optimized crossing point shifts to smaller values of p. 
However, note that the above analysis does not account 
for systematic deviations from scaling. Similar to conven- 
tional data-collapse techniques it analyses the data "as 



observed" , implying that the scaling assumption holds. 



Appendix B: Results for disk-shaped cluster growth 

Assuming a circular shape of the growing clusters (see 
Fig- EJ), and upon analysis of the order parameter (again 
for square systems with side- length up to L = 1024), we 
yield the following results for the critical properties of 
the 2D "Touch and stop" cluster growth model: 

e. Binder ratio: Considering the Binder ratio, the 
best data collapse (obtained in the range e G [—2,2]) 
yields p c = 0.4978(3), and v = 1.32(4) with a quality 
S = 0.984. 

/. Order parameter: Considering the relative size of 
the largest cluster found in the final configuration, the 
best data collapse (obtained in the range e G [—1,1]) 
yields p c = 0.4974(1), v = 1.29(3), and (i = 0.145(4) 
with a quality S = 0.29. 

g. Order parameter fluctuations: A best data col- 
lapse for the order parameter fluctuations (attained in 
the range e G [—1.5,1.5]) results in the estimates p c = 
0.4978(5), v = 1.33(2), and 7 = 2.38(3) with a quality 
S = 0.89. 

h. Cluster size distribution: At the critical point, 
the algebraic decay of the probability mass function 
n s for clusters of size s is governed by the exponent 
t = 2.050(6). This estimate was obtained for the data 
recorded at L = 1024, using a least-squares fit to the 
interval s G [500, 7000] (as before, the errorbar was esti- 
mated via bootstrap resampling). 

The observation that the critical points for rhombic 
and disk-shaped cluster growth differ only slightly can 
be understood from the time evolution of the individual 
clusters. At short times (i.e. t < 2) both "growth styles" 
feature identical clusters. Only at times t > 2, the disk- 
shaped growth style leads to more convex clusters that 
assume a circular shape in the limit t — > 00. Close to the 
critical points of both growth styles, the cluster growth 
procedure stops after maximally 4 time steps (i.e. at t = 
4). Thereby, the majority of initial seeds grow for less 
than 4 time steps and become inactive in a state where 
one cannot tell apart clusters that were grown using a 
rhombic or disk-shaped growth rule. Note that this is 
due to the discreteness of the underlying lattice. 
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